A cosmological constant from the QCD trace anomaly? 
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According to recent astrophysical observations the large scale mean pressure of our present universe 
is negative suggesting a positive cosmological constant like term. This article addresses the question 
of whether non-perturbative effects of self-interacting quantum fields in curved space-times may yield 
a significant contribution. Focusing on the trace anomaly of quantum chromo- dynamics (QCD), a 
preliminary estimate of the expected order of magnitude yields a remarkable coincidence with the 
empirical data, indicating the potential relevance of this effect. 
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Recent measurements of the cosmic microwave back- 
ground [l] suggest that the large scale structure of our 
universe is quite accurately described by the conformally 
flat Friedmann-Robertson- Walker (FRW) metric 



ds 2 = dr 2 - n 2 (r) dr 2 = fl 2 (t) (dt 2 - dr 2 



(1) 



with r denoting the proper and t the conformal time, 
respectively. The temporal variation of the scale factor 
Q 2 inducing the cosmological red-shift is represented by 
the Hubble parameter 
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Inserting the FRW metric in Eq. ([!]) into the Einstein 
equations (with the cosmological constant A) 



(3) 



the aforementioned observations and supernova J2[] data 
consistently yield the following conclusions: the 00- 
component of the r.h.s. of Eq. (^) equals (at least approx- 
imately) the critical density g ~ g cr n and the spatial ii- 
components - associated with the pressure p - are nega- 
tive: p/gcrit ~ —2/3. As a result the universe is presently 
undergoing an accelerated expansion as approximately 
described by the de Sitter metric fl(r) — exp{$jr}. 

A negative pressure together with a positive energy 
density necessarily implies a non-vanishing trace of the 
energy-momentum tensor T^. Since T M „ can be derived 
via the variation of the action A with respect to the met- 
ric Q^, i.e. T M „ = 2(-q)- 1 / 2 5A/5q^, its trace corre- 
sponds to the change of A under the conformal transfor- 
mations Q^ v (x) — > CI 2 {x)Qfiv (x) by virtue of Euler's law: 
TP = -Qi-gy^SA/SQ, see e.g. §. 

Let us focus on the contribution of the SU (3)-color 
gauge field theory of QCD to the r.h.s. of the Einstein 
equations (|j|) in the following. Its dynamics are governed 
by the well-known Lagrangian density 



£ = -- G^Gr + $ (*7^ +gT a >fAl - m) vb , (4) 



with G£„ = d^Al ~ d v Al+gf\ c A\Al being the gluonic 
field strength tensor. Here f a bc denote the SU (3) struc- 
ture constants, T a its fundamental generators, and g is 
the strong coupling. For simplicity we drop the ghost 
fields as well as the gauge fixing terms and consider only 
one single flavor, i.e. quark-species ip. The remaining 
electro-weak sector of the standard model will be dis- 
cussed at the end of this article. 

On the classical level all gauge field theories as de- 
scribed by Eq. (Q) are conformally invariant (for m = 0). 
According to the above arguments this feature implies a 
vanishing trace of the classical energy-momentum tensor 
(as one can easily check by an explicit calculation.) Turn- 
ing to the quantum field theoretical description the situ- 
ation becomes more complicated. In the first place, the 
naive expectation value of the operator-valued energy- 
momentum tensor diverges due to the infinite zero-point 
energy. In order to renormalize this singularity by an 
appropriate counter-term one has to interpret the cos- 
mological constant A in Eq. (||) as a bare quantity j|. 

After such a minimal subtraction procedure the trace 
of the renormalized expectation value of the energy- 
momentum tensor (Tj^rcn vanishes (again assuming 
m = 0) for free (g = 0) fields in flat (fl = 1) space- 
times - but not in the general case. This phenomenon is 
called the trace anomaly and goes along with the dynam- 
ical breaking of the conformal invariance of the classical 
theory in Eq. (Q). It has been calculated for two lim- 
iting cases: firstly for self-interacting quantum fields in 
flat (fi = 1) space-times || and secondly for free (g = 0) 
fields in curved space-times || . In the first case one ob- 
tains B 
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<G" GT)rcn + (l+7 m )m(#}re 
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The Callan-Symanzik or Gell-Mann-Low /3-function de- 
scribes the scale dependence (3 = [idg/d/j, of the renor- 
malized coupling <?(/z) and reflects the dynamical break- 
ing of the conformal invariance of the classical theory 
(dimensional transmutation). Similarly, the 7™ l -function 
corresponds to the running of the renormalized mass. 
The expectation values (G^ u G£ u ) rcn and (tptp) ren occur- 
ring in Eq. (||) represent the so-called gluonic and quark 
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condensates, respectively, see e.g. Hj. These inherently 
non-perturbative quantities again reflect the dynamical 
breaking of the classical scale invariance. Both are of 
great experimental relevance and their values have been 
confirmed within several contexts, see e.g. [pi. Since the 
symmetry breaking scale Aqcd is (for m = 0) the only 
scale in the theory (@), it yields (G£„GT)rcn - 0(A QCD ) 
and (^V)rcn = C(A-q CD ) - at least for m -C Aqcd- As 
it is well-known, the /3-function occurring in Eq. (|5|) can 
be calculated within the framework of perturbation the- 
ory and it turns out to be negative. Consequently it 
is now commonly accepted (cf. ||) that the QCD trace 
anomaly gives rise to a negative energy density (since 
(G^G^) rcn > and (tpi>) len < 0) of the QCD vacuum 
in the Minkowski space-time. 

However, such a huge amount of negative energy den- 
sity of order C(Aq CD ) blatantly contravenes our obser- 
vations. This drastic and global violation of the (weak 
and dominant) energy conditions (see e.g. H) in the 
Minkowski space-time goes along with a fundamental 
contradiction if one includes gravity since the r.h.s. of 
the Einstein equations (||) associated with a flat space- 
time vanishes. 

Consequently, regarding the Einstein equations (j|), 
one is led to absorb the aforementioned energy density by 
renormalizing the cosmological constant A, in complete 
analogy with the case of the zero-point energy (which 
determines the divergent part of A only). In the same 
manner as one adjusts the mass-counter term in the self- 
energy renormalization of the electron, for example, one 
has to fix the bare cosmological constant by demanding 
that the r.h.s. of the Einstein equations (^) vanishes for 
the Minkowski vacuum 



Minkowski vacuum 



0. 



(6) 



On the other hand, a non-trivial geometry of the space- 
time may also induce a non- vanishing trace - even for free 
fields (second limiting case g = 0). In this case (Tp rcn 
is given by the sum of a bilinear form of the curvature 
tensor (such as D^W" or 5H 2 ) and second derivatives of 
it [am], cf. §]. For the free QCD field (with g = m = 
0) within the de Sitter space-time (1(t) = exp{f)T}, for 
example, one finds (h = c = 1) 



(TP 
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In contrast to the contribution in Eq. (^J), there is no 
reason to absorb this term by renormalization of A. In 
view of its potentially space-time dependent character 
such a procedure would be rather strange. However, here 
the associated energy density is far too small to explain 
the observations 

In summary, the (renormalized) expectation value of 
the energy-momentum tensor acquires an anomalous 
trace for self-interacting quantum fields in flat space- 
times (||) on the one hand as well as for free fields in 



curved space-times (0) on the other hand. However, both 
effects taken alone are not capable of explaining the neg- 
ative pressure as suggested by the observations |l],|j . But 
this is just what one might expect, since realistic inves- 
tigations have to involve both contributions simultane- 
ously, i.e. the (non-perturbative) effects of self-interacting 
fields in curved space-times. A rigorous derivation of the 
renormalized expectation value of the energy-momentum 
tensor for this scenario appears to be rather involved and 
is not the aim of the present article. Here we just give a 
preliminary estimate of the expected order of magnitude 
of the effect. To this end we employ an adiabatic approx- 
imation (cf. ||) by exploiting the huge difference of the 
involved time scales. The cosmic evolution - governed 
by io - is extremely slow compared to the typical fluc- 
tuations of the quantum field as determined by Aqcd- 
Consequently the adiabatic approximation is an expan- 
sion in the small parameter Sj/ Aqcd = O(10~ 40 ). So the 
zcroth-order term is the pure Minkowski (flat space-time) 
contribution whereas the first-order term represents the 
lowest correction induced by the cosmic expansion. 

In order to calculate the renormalized expectation 
value of the energy-momentum tensor it is essential to 
specify the correct vacuum state associated with our 
expanding universe, cf. S. To this end we adopt the 
Schrodinger picture 



dt 



I*) = -^FHW(0I*). 



(8) 



where -ffFRwW denotes the Hamilton operator, i.e. the 
generator of the time evolution, with respect to the con- 
formal coordinates (t, r) in Eq. ([!]). Within the adiabatic 
approximation, the explicitly time-dependent Hamilto- 
nian iJpRw(i) of an expanding universe can be related 
to the (time-independent) Minkowski Hamiltonian -ffMin 
via 

HFnw(t) = exp j— iO^iSj # M iiiexp {+i(l(t)s} , (9) 

with S being the generator for the conformal transfor- 
mations g M „ — > d 2 Q^ v in the Schrodinger picture. In 
terms of a dynamically scaled state defined via \iff) = 
exp{+i(}(t)S}\^>} , the Schrodinger equation assumes the 
form 



(10) 



-\$) = -i[H mil -ns)m 



Treating Hi = —QS = 0(Sj) as a perturbation 
and switching to the interaction representation S(t) = 

expj+i/fMini} S exp{— iHuint} we may solve the above 
equation in linear response, i.e. first order adiabatic ex- 
pansion 

o 

|$) - |$ in ) + i f dt (1(f) S(t) |* in ) + O(rf) . (11) 
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If we assume f2(£) = £lo + cxp{Sjt}, it is reasonable to 
take the Minkowski vacuum |0Min) with -HiviinlOMin) = 
as the initial condition I'I'in) = Omui)- The remaining 
time integration yields H^ in and hence we arrive at 



|0 F RW) = |0 M i„) + Si H^ln S |0 M in) + 0(Si 



(12) 



Therefore the adiabatic QCD vacuum |0frw) of an ex- 
panding universe is not the instantaneous ground-state 
|0Min) of i?FRw(£) or Hmiu ~ it acquires corrections al- 
ready in the first order of the adiabatic expansion. In- 
stead it is the ground-state of the corrected Hamiltonian 
Hwrn — QS, which can easily be verified using stationary 
perturbation theory (see e.g. || for free fields). 

The remaining question is, of course, whether the first- 
order correction fj H^l n S OMin} to the vacuum state en- 
tails a first-order correction to the expectation value of 
T^ v . In order to illustrate this point let us consider the 
simple example of a time-dependent harmonic oscillator 



H(t) 



2fi(t) 



p 2 + n 2 (t)Q : 



(13) 



where uj corresponds to Aqcd ■ In examining the question 
of whether the dynamical scale symmetry breakdown in 
QCD can be modeled by such a simple quadratic po- 
tential one might consider the CP(N — 1) or the O(N) 
cr-models |^|. These strongly interacting theories repro- 
duce several features of QCD, such as dynamical scale 
symmetry breakdown. They can be solved in the large 
iV-limit and in the leading order they effectively behave 
like massive free fields. After a normal mode decom- 
position one therefore indeed obtains terms like the one 
above. For the Hamiltonian in Eq. (|TJ) a change of the 
scale factor Q(t) as in Eq. (||) can simply be generated 
by the squeezing operator 



(14) 



Consequently the expectation values of operators such as 
Q2 Qr p2 ^ n0 £ aC q U j re a first-order correction. This re- 
sult can be transfered directly to free quantum fields: 
pitching on a particular normal mode with the wave- 
length k the conformal charge S again acts like a squeez- 
ing operator S — » i(a' k a_ k — cikci-k)- As a result there is 
no first-order correction to the expectation value of T^„ 
for free (linear) fields. 

However, if we leave the free-field sector and take in- 
teractions into account the situation may change: let us 
consider the interaction Hamiltonian 



H int (t)=g f dt'G rct (t-t'){Q 2 (t),Q 2 (t')} , (15) 



where Q(i) — Q cos(ujt) + P sm(ojt) / w denotes the unper- 
turbed time-dependent operator in the interaction pic- 
ture. The retarded propagator G rct (t — t') encodes the 



dynamics of an intermediate (interaction) degree of free- 
dom which has been integrated out. Again such a term 
can be motivated by the CP(N — l)-models: in the large 
TV-limit these strongly interacting massless theories ef- 
fectively transform into massive fields obeying weak (i.e. 
next-to-leading order in l/N) long-range four-point in- 
teractions, cf. ||. 

For general Green functions G xe t(t — £'), the operator 
in Eq. ( |l5| ) does entail a first-order correction. This can 
be most easily verified by assuming g <c 1 which allows 
for a perturbative treatment. 

In view of these considerations one might expect a first- 
order contribution to (T M „) rcn to be possible in the case 
of QCD: since the classical as well as the free quantum 
field in Eq. (^) are (for m = 0) conformally invariant, 
their solutions would simply be scaled during the expan- 
sion of the universe - like the red-shift of the photon 
field. (This would actually happen if the universe were 
to expand very rapidly Sj ^> Aqcd-) However, the strong 
self-interaction on the quantum level breaks the confor- 
mal invariance and introduces a fixed scale Aqcd lead- 
ing to a positive pressure given by Eq. (|5|). Within an 
expanding universe the balance of these two tendencies, 
i.e. following the expansion on the one hand and retaining 
the scale on the other hand, leads to a displaced vacuum 
state dl|). 

Let us assume that a part of the positive vacuum pres- 
sure in Eq. (|^) can be explained by relatively localized 
(non-perturbative) vacuum fluctuations (e.g. instantons 
or oscillons ||) which repel each other (at least in 
average, cf. 0). Let us further assume that the dynami- 
cal breakdown of the scale symmetry is basically encoded 
by these (non-perturbative) vacuum-fluctuations whereas 
their (repulsive) interactions are adequately described by 
the free (perturbative) and thus conformally invariant 
field equations. In this case their solutions would sim- 
ply be scaled during the cosmic expansion in contrast 
to the non-perturbative fluctuations which retain their 
scale and hence are not affected. Within an expanding 
universe, then, every vacuum fluctuation "sees" all other 
vacuum-fluctuations "red-shifted" , i.e. their repulsion ac- 
quires a correction proportional to SiR, where R denotes 
their (mean) distance (cf. 0). Accordingly, the positive 
vacuum pressure in Eq. (||) gets diminished by an amount 
of first order in Si. 

Based on this intuitive picture it appears plausible to 
admit a correction to the expectation value of the energy- 
momentum tensor within the FRW vacuum in Eq. (|12| ) 
which is linear in Sj. After the renormalization described 
in Eq. i.e. the subtraction of the Minkowski contri- 
bution, we therefore obtain 



Irfi \FRW 



— Si ( H Min 
= 0(£A 3 Q cd) 



S) 



Min 

ren 



H.c. + 0(Si 2 



(16) 



Let us estimate the associated order of magnitude: al- 
though Aqcd depends on the renormalization scheme 
we may fix it approximately via Aqcd = O(10 8 eV) = 



3 



O(10 14 m^ 1 ). The masses of the light quarks which 
dominantly couple to the gluonic field are roughly of 
a similar order of magnitude. The Hubble expansion 
parameter $j is about 10 _26 m _1 . Inserting the above 
values we finally arrive at (T M „)™ W = 0(lO 16 m- 4 ) or 

(2>)™ W = 0(ICT 29 gcm- 3 ). By inspection one finds 
that the deduced order of magnitude nicely fits the em- 
pirical data £> cr it ~ 10~ 29 gcm~ 3 . In view of the huge 
difference of the involved scales [5}/Aqcd = O(10 -40 )] 
this remarkable coincidence seems to be almost too good 
to be just an accident. At least it indicates the poten- 
tial relevance of the effect described in the present arti- 
cle with regard to the interpretation of the astrophysical 
data 

It should be mentioned here that a pressure induced by 
the expansion of our universe with p oc fj generates a cos- 
mic evolution which differs from that with a true cosmo- 
logical constant p = const: by inserting the FRW metric 
(|l|) into the Einstein equations (||) one obtains the Fried- 
mann equation 3Sj 2 = 8ttG^q. Furthermore the Einstein 
equations imply V /J (T /ily } rcn = 0, i.e. dg/dt = —3(g+p)$j. 
Combining these two equalities and specifying the pres- 
sure p one may determine the time-evolution of our uni- 
verse. Unfortunately the presently available data (such 
as the Hubble parameter or the age of the universe) are 
not precise enough to distinguish the two cases (p oc Sj 
and p = const) . 

Of course one may ask whether the remaining electro- 
weak sector of the standard model generates similar con- 
tributions: typically (see e.g. j^]) non-perturbative ef- 
fects (such as (G" G% u ) len and (V^ren) display a depen- 
dence on the coupling of exp{— 8n 2 /g 2 } = exp{— 2n/a}. 
The scale of the dynamical symmetry breaking Aqcd 
obeys a similar non-analytical dependence on the cou- 
pling g. Inserting o;qed ~ 1/137 into the above expres- 
sion one obtains a suppression by an order of magnitude 
of 10~ 370 . Hence the contributions arising from the dy- 
namical breaking of scale invariance can safely be ne- 
glected in this case. 

The remaining explicit breaking of the scale symme- 
try induced by the Higgs field of course also generates 
contributions to T^ v . However, the general structure of 
all these terms is given by m 2 ($^$} ron and according 
to the arguments after Eq. ( |l4[) they do not contribute 
to the first order in f). Although the mixture of these 
terms caused by interactions remains subject to further 
considerations, a contribution of the electro-weak sector 
in analogy to QCD is not obvious. 

In summary, the present article motivates a deeper ex- 
amination of the vacuum of strongly interacting fields in 
the gravitational background of our expanding universe 
- for the present epoch as well as for earlier stages, cf. || . 
These investigations might perhaps lead to a better un- 
derstanding of some of the problems in cosmology with- 
out necessarily invoking yet unknown low-energy fields, 
for example quintessence (see e.g. [To|]). 
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